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A  TWO-STROKE  DIESEL  ENGINE  SIMULATION  PROGRAM 

Jon  H.  Van  Gerpen 
Engine  Research  Institute 
Iowa  State  University 
Ames,  Iowa 


SUMMARY 

A  computer  program  simulating  a  two-stroke  diesel  engine  is  developed  and 
documented.  The  program  is  suitable  for  simulating  the  diesel  core  of  a  high- 
output  combined-cycle  diesel  engine.  The  engine  cylinder  and  the  intake  and 
exhaust  ports  are  defined  as  independent  thermodynamic  systems  and  the  mass  and 
energy  equations  for  these  systems  are  developed.  A  single  zone  combustion 
model  is  used  and  perfect  mixing  during  scavenging  is  assumed.  The  program 
input  requirements  and  output  results  are  discussed.  A  sample  case  is  provided 
for  an  opposed  piston,  uniflow  scavenged  two-stroke  diese1  engine.  7 r“/,  '■ 


INTRODUCTION 


This  report  documents  a  computer  program  solving  the  equations^from  a 
mathematical  model  that  simulates  a  two-stroke  diesel  engine.  The  model  con¬ 
siders  the  thermodynamics  and  fluid  mechanics  of  the  working  fluid  from  the 
entrance  of  the  intake  port  to  the  exit  of  the  exhaust  port.  The  program  car/ 
predict  the  effect  on  engine  performance  of  changes  in  parameters  such  as 
speed,  boost  pressure,  valve  timing,  and  fueling  level.  The  program  output 
provides  information  about  power  output,  brake  mean  effective  pressure  (BMEP) , 
heat  transfer  losses,  and  cylinder  pressures  and  temperatures.  The  program  was 
written  in  modular  form  so  that  the  submodels  could  be  modified  or  replaced 
without  requiring  program  alteration. \ 


/ 

Most  of  the  Input  data  is  supplied  by  an  input  data  set,  but  because  the 
program  is  Intended  to  be  used  as  part  of  a  larger  program  that  will  match  the 
diesel  with  attached  turbomaxW nery,  the  inlet  and  exit  conditions  and  the  mass 
f low a  subroutine  argument  list. 

‘-This  report  discusses  the  mathematical  model  used  in  the  program  and 
describes  the  subroutines  that  make  up  the  model.  Input  requirements  are 
stated  and  the  output  listing  is  explained.  A  sample  program  run  is  also 


provided,  f 


SYMBOLS 


ATDC  after  top  dead  center 


C1  shape  parameter  for  diffusion  burning 

C2  crank  angle  for  start  of  combustion 

combustion  duration  parameter 

STATEMENT  "A"  per  Robert  Bill 
Propulsion  Pirrctoratc/AVSCOM,  US  Army 
Aviation  Research  and  Tecnology  Activity 
Cleveland,  OH  44135 
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FC 

f 


f  ,tot 
fuel  ,new 
IP 

ri 

[2 
mR,BA 

mR,fuel 

mR,UA 

P 

P 

• 

0 


premixed  burning  fraction 
shape  factor  for  premixed  burning 
stoichiometric  fuel-to-alr  ratio 
residual  fraction 
enthalpy 

enthalpy  of  gases  exiting  to  exhaust  manifold 
enthalpy  of  gases  exiting  to  exhaust  port 
enthalpy  of  gases  entering  intake  port 
enthalpy  of  gases  entering  cylinder 
mass 

mass  of  cyl inder  gases 

mass  flow  rate  of  gases  from  exhaust  port  to  exhaust  manifold 

mass  flow  rate  of  gases  from  cylinder  to  exhaust  port 

mass  of  fresh  charge  (pure  air) 

mass  flow  rate  of  fuel  added  to  cylinder 

total  amount  of  fuel  added  to  cylinder 

mass  of  new  fuel  Injected 

mass  at  Intake  port 

mass  flow  rate  of  gases  from  Intake  manifold  to  Intake  port 

mass  flow  rate  of  gases  from  Intake  port  to  cylinder 

mass  of  burned  air  In  residual 

mass  of  burned  fuel  In  residual 

mass  of  unburned  air  In  residual 

purl ty 

pressure 

heat  transfer  rate 
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R  gas  constant 

Rjp  gas  constant  at  Intake  port 

rQ  del i very  ratio 

s  entropy 

T  temperature 

Tjp  temperature  at  intake  port 

1)  ^  internal  energy  of  the  cylinder  gases 

UEp  internal  energy  of  the  exhaust  port  gases 

Ujp  internal  energy  of  the  intake  port  gases 

u  internal  energy  per  unit  mass 

uc^  internal  energy  of  the  cylinder  gases  per  unit  mass 

V  volume 

W  rate  of  doing  work,  also  equal  to  p  dV/de 

y  nondimenslonal  time  parameter 

<j>cyl  equivalence  ratio  of  cylinder  gases 

4>c  equivalence  ratio  of  mass  leaving  exhaust  port  to  exhaust  manifold 

tl 

4>c  equivalence  ratio  of  mass  leaving  cylinder  to  exhaust  port 

t2 

4>t  equivalence  ratio  of  mass  entering  cylinder  from  Intake  port 

2 

e  crank  position,  or  crank  angle 

nch  charging  efficiency 

nsc  scavenging  efficiency 

ntr  trapping  efficiency 

THE  MATHEMATICAL  MODEL 

The  mathematical  model  for  the  engine  is  based  on  applying  the  first  law 
of  thermodynamics  to  the  systems  composed  of  the  engine  cylinder  and  the  Intake 
and  exhaust  ports.  The  ideal  gas  equation  and  equilibrium  gas  properties  are 
also  used  to  solve  for  the  pressure  and  temperature  in  the  cylinder. 
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Figure  1  shows  the  systems  being  analyzed  in  this  section.  The  first  sys¬ 
tem,  the  engine  cylinder,  receives  mass  from  the  Intake  port  system  and  deliv¬ 
ers  mass  to  the  exhaust  port  system.  The  possibility  for  backflow  exists  when 
the  cylinder  pressure  is  above  the  Intake  port  pressure  or  below  the  exhaust 
port  pressure.  The  first  law  of  thermodynamics  can  be  written  as  follows  for 
the  system  consisting  of  the  engine  cylinder: 


dU 


cyl 


de 


=  Q  -  W  + 


I  "I 
2  l2 


m^h^ 


m 


E  E 
t2  2 


(1) 


where 


UCyi  the  Internal  energy  of  the  cylinder  gases 

Q  the  heat  transfer  rate  from  the  cylinder 

• 

W  the  rate  at  which  work  is  done  by  the  cylinder  gases 

mT  the  mass  flow  rate  Into  the  cylinder  from  the  Intake  port 

2 

mr  the  mass  flow  rate  from  the  cylinder  to  the  exhaust  port 
fc2 

mf  the  mass  flow  rate  of  fuel  added  to  the  cylinder 

hT  enthalpy  of  gases  entering  from  Intake  port 

2 

hf  enthalpy  of  fuel  added  to  the  cylinder 

hF  enthalpy  of  gases  exiting  cylinder  to  exhaust  port 

t2 

9  crank  position,  or  crank  angle  (used  as  a  time  parameter) 
The  derivative  of  Uc^  can  be  written  as 


dU 


cyl 


d0 


d<"1cVILlcvl> 

de 


d<ucvl> 


dm 


+  u. 


cyl 


(2) 


"'cyl  d9  T  ucyl  de 

where  nvy]  is  the  mass  of  cylinder  gases  and  uCyi  is  the  specific  Internal 
energy  or  the  cylinder  gases. 

The  Ideal  gas  equation  can  be  differentiated  with  respect  to  crank  angle 
to  give  the  following  expression: 


dV  \i  mD  dT  mT  dR  DT  dm 
+  V:rs-mR-7^  +  mT-^r+RT-7^ 


de 


de 


de 


de 


de 


(3) 


Each  of  the  properties  p,  V,  m,  R,  and  T  refer  to  the  cylinder  gasrt: 
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Derivatives  of  u  and  R  with  respect  to  9  can  be  expressed  In  terms 
of  derivatives  of  T,  P,  and  <J>,  the  equivalence  ratio,  as  shown  In  equations 
(4)  and  (5).  The  equivalence  ratio  used  here  Is  defined  as  the  actual  fuel-to- 
alr  mass  ratio  divided  by  the  stoichiometric  fuel-to-alr  ratio. 

dJR  3R  d]3  3R  dT  3R  d£_ 

de  =  ap  de  +  3T  de  +  a<t>  de 

du  3u  d£  9u  dT  3u  d£_ 

de  =  ap  de  aT  de  +  a$  de 


Equations  (3)  and  (4)  can  be  combined  and  solved  for  dT/d6.  This  equation  can 
then  be  used  to  eliminate  dT/de  In  the  equation  that  results  from  combining 
equations  (1)  and  (5)  to  obtain  an  expression  for  dp/de. 

The  equivalence  ratio  In  the  cylinder  changes  because  of  the  addition  of 
air  from  the  Intake  port,  fuel  from  the  fuel  Injector,  or  backflow  from  the 
exhaust  port.  The  following  expression  can  be  obtained  for  the  variation  of 
4>Cy1  with  crank  angle: 


+  ^cylFs 


*12  '  *cyl  , 


♦e2  -  *cyl 

'  +  *C  Fc 


(6) 


where  *cy) 


are  defined  previously  and 


where 


<j>j  equivalence  ratio  of  mass  entering  or  leaving  cylinder  through  the  Intake 
l2  opening 

<|>E  equivalence  ratio  of  mass  entering  or  leaving  cylinder  through  the 
l2  exhaust  opening 

F$  stoichiometric  fuel-to-alr  ratio 

The  equations  described  previously  for  dp/de,  dT/d0,  and  d4»/d0  can  be 
Integrated  simultaneously  to  obtain  the  p,  T,  and  $  In  the  cylinder  as  a 
function  of  0.  The  heat  transfer  rate  Is  calculated  with  Annand's  correlation 
<r*f.  1).  The  mass  entering  and  leaving  the  cylinder  through  the  Intake  and 
exhaust  openings  Is  calculated  on  the  basis  of  an  Isentroplc  expansion  Into  a 
region  of  lower  pressure.  Deviations  from  this  Ideal  flow  are  Incorporated 
through  the  use  of  a  user-supplied  flow  coefficient.  The  mass  flow  rate  of 
fuel  added  to  the  cylinder  Is  also  user-specified  through  a  functional  expres¬ 
sion  given  In  subroutine  DWCA. 


5 


This  model  assumes  that  the  mass  leaving  a  system  has  the  same  properties 
as  the  system.  This  assumption  implies  that  there  will  be  perfect  and  instan¬ 
taneous  mixing  between  any  mass  added  to  the  system  and  the  mass  already  exist¬ 
ing  in  the  system.  This  has  implications  for  the  scavenging  process  because 
the  actual  scavenging  process  can  be  more  or  less  effective  than  perfect  mix¬ 
ing.  If  there  is  considerable  short  circuiting  of  the  intake  air  directly  to 
the  exhaust  port,  then  the  actual  scavenging  efficiency  will  be  less  than  cal¬ 
culated  here,  but,  if  the  incoming  air  is  able  to  push  the  exhaust  products  out 
of  the  cylinder  without  excessive  mixing,  then  the  actual  scavenging  efficiency 
will  be  higher.  The  actual  scavenging  efficiency  in  an  engine  is  highly  sensi¬ 
tive  to  combustion  chamber  design  and  engine  operating  conditions,  so  the  per¬ 
fect  mixing  assumption  represents  a  reasonable  compromise  that  should  still 
allow  the  prediction  of  trends. 

The  second  two  systems  identified  in  figure  1,  the  intake  port  and  the 
exhaust  port  are  described  by  equations  similar  to  those  developed  previously. 
The  following  assumptions  apply  to  these  systems: 

(1)  Volumes  of  intake  and  exhaust  ports  are  constant. 

(2)  Pressures  in  intake  and  exhaust  ports  are  constant  and  equal  to  the 
user-specified  values  of  intake  and  exhaust  manifold  pressures. 

(3)  The  intake  and  exhaust  ports  are  adiabatic.  There  is  no  heat 
transfer. 

(4)  Mass  passing  from  the  intake  manifold  into  the  intake  port  is  assumed 
to  be  air  only.  No  combustion  product  residual  gases  are  present.  If  backflow 
occurs  from  the  intake  port  system  to  the  intake  manifold,  this  mass  is  also 
assumed  to  be  air  only,  regardless  of  the  actual  equivalence  ratio  in  the 
Intake  port. 

(5)  Mass  passing  from  the  exhaust  port  to  the  exhaust  manifold  has  proper¬ 
ties  equal  to  the  current  values  In  the  exhaust  port.  If  backflow  occurs  from 
the  exhaust  manifold  to  the  exhaust  port,  this  mass  also  has  properties  equal 
to  the  current  properties  of  the  gases  in  the  exhaust  port. 

The  first  law  of  thermodynamics  can  be  written  as  follows  for  the  Intake 
port  system: 


n,  hr  —  m*  hj 
L]  X1  LZ  x2 


(7) 


where 

Ujp  the  internal  energy  of  the  intake  port  gases 

m,  the  mass  flow  rate  into  the  Intake  port  from  the  Intake  manifold 
X1 

tu  the  enthalpy  of  the  gases  entering  the  Intake  port  from  the  intake 
*1  manifold 
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The  next  step  is  to  differentiate  the  ideal  gas  equation  subject  to  the 
constraint  that  the  pressure  and  volume  are  constant. 


dT 


0  =  mIPRIP 


IP 


d6 


dRjp  dm^p 

+  mIPTIP  ~W  +  RIPTIP  ~W 


(8) 


Equations  (4)  and  (5)  can  be  used  to  represent  the  derivatives  of  u  and  R 
with  respect  to  9  although  the  pressure  terms  can  be  eliminated  because  of 
the  assumption  that  the  intake  port  pressure  is  constant. 

The  equivalence  ratio  in  the  intake  port  changes  both  because  of  the  addi¬ 
tion  of  air  from  the  intake  manifold  and  because  of  a  possible  backflow  from 
the  cylinder.  The  following  expression  can  be  obtained  for  the  variation  of 
<J>IP  with  crank  angle: 


d4 ’IP 

(I  +  *IPFS) 

'(♦ip  ■  *I2) 

de  = 

mIP 

1  *  *I2Fs 

<t»ipmr 


(9) 


where 

<j)jp  equivalence  ratio  of  mass  in  intake  port 

In  a  manner  similar  to  the  case  for  the  cylinder,  equation  (8)  can  be 
solved  for  dT/d9  to  obtain  the  differential  equation  for  temperature;  this 
equation  can  then  be  used  to  eliminate  the  dT/d©  term  that  appears  in  equa¬ 
tion  (7)  after  substitution  of  equation  (5).  In  this  case,  the  unknown  in 

equation  (7)  becomes  m^  .  Thus,  the  equations  for  mIp,  dTIp/d0,  and  d<j>Ip/de 

can  be  integrated  along  with  the  equations  for  the  in-cylinder  conditions  to 
obtain  mjp,  Tjp,  and  4>i p .  The  intake  port  conditions  are  of  little  interest 
themselves,  but  they  are  required  to  fix  the  boundary  conditions  of  the  equa¬ 
tions  describing  in-cylinder  conditions.  Although  they  will  not  be  listed 
here,  equations  can  be  developed  for  the  exhaust  port  in  exactly  the  same  man¬ 
ner  as  for  the  intake  port. 

The  equations  described  previously  must  be  integrated  through  the  engine 
cycle.  Additional  quantities,  while  not  required  to  simulate  the  system,  will 
also  be  summed  during  the  cycle.  There  are  a  total  of  17  equations.  The  inte¬ 
gration  process  is  started  at  the  crank  position  where  the  exhaust  valve  opens. 
Initial  values  are  required  for  each  of  the  variables.  While  most  of  the  vari¬ 
ables  are  initialized  to  zero,  variables  such  as  the  cylinder  temperature  and 
pressure  must  have  initial,  nonzero  values.  Since  these  values  will  not  gener¬ 
ally  be  known  in  advance,  values  are  assumed  and  the  integration  is  started. 
Since  the  process  being  described  is  a  mechanical  cycle,  the  properties  at  the 
end  of  the  integration  should  be  the  same  as  at  the  start.  The  solution  pro¬ 
cess  is  Iterative,  where  the  final  values  of  the  Integrated  variables  are  used 
as  the  initial  values  of  the  next  Iteration.  When  the  difference  between  the 
initial  and  final  values  is  less  than  a  prescribed  value,  the  program  is  con¬ 
sidered  to  have  converged  and  the  program  halts. 
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The  program  offers  two  options  for  integration  routines.  A  subroutine 
called  RKGS  can  be  used  and  is  supplied  with  the  DSL2  program.  However,  if 
the  IMSL  library  is  available,  the  fifth-order  Runge-Kutta  integration  routine 
DVERK,  which  can  be  used  instead,  results  in  much  faster  execution  time. 
Statements  to  call  each  routine  are  provided  in  the  program. 


SUBROUTINE  DESCRIPTIONS 


INPUT1:  This  subroutine  reads  in  the  data  set  containing  the  operating  param¬ 
eters  for  the  engine  being  simulated.  A  description  of  the  input  quantities  is 
contained  in  the  section  on  input  data.  This  subroutine  also  does  some  prelim¬ 
inary  processing  of  the  input  data.  In  particular,  this  subroutine  calls  the 
SPLINE  subroutine  to  generate  the  array  of  spline  coefficients  for  the  intake 
and  exhaust  flow  areas.  The  flow  area  information  is  supplied  with  the  input 
data,  either  as  a  table  of  areas  versus  crank  position  or  as  port  geometry,  and 
then  a  cubic  spline  is  fit  to  this  data.  This  operation  is  only  performed  once 
and  after  this  time,  the  spline  coefficient  array  is  accessed  by  the  function 
subroutines  AREAIN  and  AREAEX  to  calculate  the  intake  and  exhaust  flow  areas 
for  any  crank  angle. 


MODEL:  This  subroutine  evaluates  the  derivatives  of  the  differential  equations 
that  result  from  the  mathematical  model  of  the  diesel  engine.  MODEL  determines 
whether  the  intake  and  exhaust  valves  are  open,  and  if  so,  calls  TMI2  and  TME2 
to  get  the  mass  flow  rates  into  and  out  of  the  cylinder.  Then  it  calls  PERX  to 
get  the  properties  of  the  gases  entering  and  leaving  the  cylinder.  These  quan- 

setting  up  the  derivatives  of  the  first-order  differential 
be  integrated  by  DVERK.  The 
the  subroutine  in  the  Y  array 
tion  derivatives  are  returned  in  the  DERY  array. 

17  simultaneous  equations.  The  quantities  being 
follows : 


tities  are  used  in 
equations  that  wi 1 1 
tion  values  enters 


latest  estimate  of  the  func- 
and  the  newly  evaluated  func- 
The  program  is  integrating 
integrated  are  defined  as 


Y( 1 )  cylinder  pressure 

Y(2)  cylinder  temperature 

Y(3)  equivalence  ratio  of  cylinder  contents 

Y<4)  mass  of  gases  in  cylinder 
Y(5)  cumulative  mass  entering  cylinder 
Y(6)  temperature  of  gases  in  intake  port 
Y<7)  equivalence  ratio  of  gases  In  intake  port 
Y(8)  cumulative  mass  leaving  cylinder 
Y<9)  temperature  of  gases  in  exhaust  port 
Y< 1 0)  equivalence  ratio  of  gases  in  exhaust  port 
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Y(11)  energy  entering  cylinder  with  the  Intake  air 

Y(12)  energy  leaving  cylinder  with  the  exhaust  products 

Y ( 1 3 )  work  done  by  cylinder  gases 

Y(14)  cumulative  mass  of  fuel  injected 

Y(15)  cumulative  heat  transfer 

Y(16)  cumulative  air  entering  intake  port 

Y(17)  cumulative  air  leaving  exhaust  port 

RKGS:  This  is  a  utility  integration  routine.  It  uses  a  fourth-order  Runge- 
Kutta  technique.  RKGS  calls  the  subroutine  FCT,  which  in  turn  calls  the  sub¬ 
routine  MODEL  to  obtain  the  derivatives  of  the  differential  equations  being 
Integrated.  Although  the  program  is  currently  set  up  to  use  RKGS,  statements 
are  provided  to  change  to  DVERK,  a  routine  from  the  IMSL  library.  This  routine 
provides  much  faster  execution  times  than  RKGS. 

SPLINE:  This  is  a  utility  routine  used  to  calculate  coefficients  of  a  cubic 
spline  for  the  Intake  and  exhaust  flow  areas  versus  crank  position. 

AREAIN:  This  is  a  function  subroutine  that,  given  a  crank  angle  as  an  input 
argument,  searches  through  an  array  of  spline  coefficients  to  find  the  appro¬ 
priate  set  and  then  evaluates  the  spline  to  provide  the  Instantaneous  intake 
valve  or  port  flow  area. 

AREAEX:  This  function  subroutine  is  Identical  to  AREAIN  In  that  it  evaluates  a 
cubic  spline  to  find  the  Instantaneous  exhaust  valve  or  port  flow  area  corre¬ 
sponding  to  a  specified  crank  angle. 

PERX:  This  Is  a  utility  routine  used  to  evaluate  equilibrium  thermodynamic 
properties  of  combustion  product  gases.  Given  T,  P,  and  <{>  (the  equivalence 
ratio),  this  routine  will  compute  u,  h,  R,  s,  and  their  derivatives  with 
respect  to  T,  P,  and  <j>.  All  input  and  output  information  for  this  routine  is 
passed  through  two  COMMON  blocks.  Documentation  of  this  routine  is  provided  in 
reference  2. 

VOL:  This  subroutine  calculates  the  cylinder  volume  for  a  specified  crank 
angle.  When  the  subroutine  argument  MODE  is  set  correctly,  the  subroutine  can 
calculate  the  volume  for  a  conventional  piston-cylinder  geometry  or  for  an 
opposed  piston  configuration  with  both  crankshafts  in  phase. 

DWCA:  This  function  subroutine  is  used  to  compute  the  fuel  mass  burning  rate 
supplied  to  the  single-zone  first-law  energy  balance  for  the  cylinder  given  by 
equation  (1).  According  to  a  function  developed  by  Watson,  Pilley,  and  Marzouk 
for  simulating  diesel  combustion  (ref.  3),  the  mass  burning  rate  is  calculated. 
Although  this  function  was  developed  for  diesel  engines  operating  over  a  narrow 
range  of  conditions,  It  contains  five  parameters  that  can  be  used  to  specify 
the  shape  of  almost  any  conceivable  burning-rate  profile.  The  parameters  are 
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denoted  by  C] ,  C2,  C3,  C4,  and  C5.  The  function  describing  the  burning  rate 
Is  obtained  by  superimposing  two  algebraic  functions.  The  first  function  Is  an 
exponential  curve  given  by 


F 


1 


6.907755(C5  + 


exp 


-6.907755y 


(10) 


where 


and 


e  -  c? 

y  =  — p -  for  0  >  C 

L3 


y  =  C  for  0  <  C2 


This  function  Is  intended  to  represent  the  relatively  slow  and  controlled  dif¬ 
fusion  burning  in  the  cylinder.  A  second  curve,  given  by 


-  5000(C5 


c 

+  1  )y 


4999 


(11) 


represents  rapid,  uncontrolled  premixed  burning.  These  two  functions  are  com¬ 
bined  with  a  weighting  factor  that  determines  the  relative  amount  of  fuel  that 
burns  In  the  premixed  and  diffusion  modes 


mf  “  mf,tot[C4F2  +  (1  "  C4)Fl]  (12) 

where  mftot  Is  the  total  amount  of  fuel  added  to  the  cylinder.  Therefore, 
as  used  In  these  functions,  the  parameters  C]  and  C5  represent  shape  fac¬ 
tors  for  the  diffusion  and  premixed  burning  curves,  respectively.  C2  is  the 
crank  angle  at  which  burning  Is  Intended  to  start,  and  C3  Is  proportional  to 
the  burning  duration.  C4  Is  the  weighting  facier  that  represents  the  fraction 
of  fuel  burned  In  the  premixed  mode. 

Q:  This  function  subroutine  calculates  the  heat  loss  from  the  cylinder  gases 
to  the  walls  during  combustion.  The  calculation  follows  Annand's  correlation 
described  In  reference  1.  This  correlation  has  the  feature  that  It  treats  the 
radiation  and  convection  losses  separately.  The  cylinder  walls  are  divided 
Into  three  areas  for  the  heat  transfer  calculations:  the  head,  the  piston,  and 
the  cylinder  bore  or  sleeve.  The  head  and  piston  areas  are  constant  but  the 
sleeve  area  changes  as  the  piston  moves  up  and  down.  The  wall  temperature  data 
for  each  of  these  areas  are  entered  by  the  user  through  the  variables  THEAD, 
TPISTN,  and  TSLEEV  in  the  input  data  set.  While  the  variation  with  time  of 
the  in-cylinder  heat  transfer  Is  calculated  using  a  heat  transfer  coefficient 
determined  from  Annand's  correlation,  the  total  amount  of  heat  loss  from  the 
cylinder  will  normally  need  to  be  adjusted,  or  tuned,  by  the  user  to  give  the 


10 


expected  fraction  of  fuel  energy  lost  to  heat  transfer.  This  is  done  by  vary¬ 
ing  the  ANNND  parameter  in  the  input  data  set.  When  set  equal  to  1 .G  the  heat 
transfer  is  calculated  without  modification,  according  to  Annand's  correlation. 
If  ANNND  is  different  than  1.0,  then  it  is  used  as  a  multiplier  to  increase  or 
decrease  the  in-cylinder  heat  transfer  coefficient  that  correspondingly  changes 
the  amount  of  heat  loss. 

ISEN:  This  is  a  utility  subroutine  that  uses  PERX  to  compute  the  properties  at 
the  end  of  either  an  isentropic  compression  or  expansion.  The  calling  program 
provides  the  T,  P,  and  <j>  of  the  initial  state  and  the  final  pressure,  and 
ISEN  iteratively  determines  the  temperature  of  the  final  state  that  has  the 
same  entropy.  Although  this  routine  uses  the  assumption  of  constant  specific 
heats  to  obtain  an  initial  estimate  of  the  solution,  the  subroutine  is  intended 
*or  those  situations  where  the  specific  heats  vary  significantly  during  the 
isentropic  process. 

TMI2:  This  subroutine  calculates  the  mass  flow  rate  through  t lie  intake  valve. 
By  checking  the  pressures  in  the  intake  port  and  the  cylinder,  TMI2  determines 
whether  the  flow  is  into  or  out  of  the  cylinder.  TMI2  calls  ISEN  to  determine 
the  properties  of  the  downstream  flow  and  calls  AREAIN  to  get  the  instantaneous 
flow  area.  These  quantities  allow  the  determination  of  the  ideal  or  isentropic 
flow  through  the  valve.  Then  a  user-supplied  flow  coefficient  is  applied  to 
obtain  an  adjusted  or  actual  flow  rate  into  or  out  of  the  cylinder.  When  back- 
flow  occurs  from  the  cylinder  into  the  intake  port,  the  calculated  flow  rate  is 
given  a  minus  sign  to  denote  the  mass  loss  from  the  cylinder. 

TME2 :  This  subroutine  is  the  same  as  TMI2  except  that  the  calculations  are 
performed  for  the  exhaust  valves. 

OUTPUT:  This  subroutine  performs  some  of  the  final  calculations  and  writes 

the  results  out  both  to  the  terminal  and  to  an  output  data  file.  Friction  is 
calculated  using  an  empirically  based  polynomial  that  depends  only  on  the  mean 
pistor,  speed.  A  friction  mean  effective  pressure  is  calculated,  then  subtrac¬ 
ted  from  the  indicated  mean  effective  pressure  to  give  the  brake  mean  effective 
pressure.  This  subroutine  also  calculates  the  various  efficiencies  used  to 
characterize  the  scavenging  process.  The  calculations  are  described  in  more 
detail  in  the  section  OUTPUT. 


INPUT  DATA 

The  diesel  program  is  intended  to  be  run  as  a  subroutine  with  part  of  the 
input  variables  in  the  subroutine  argument  list.  The  subroutine  argument  list 
1  s  as  fol lows : 


SUBROUTINE  DSL2(T1 ,P1 ,P2 ,T2 ,XM) 


where 

T1 

intake 

manifold 

temperature,  °R 

INPUT 

PI 

Intake 

manifold 

pressure,  psia 

INPUT 
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INPUT 


P2  exhaust  manifold  pressure,  psia 

T2  exhaust  temperature,  °R  OUTPUT 


XM  air 

flow  rate  into  engine,  lbm/s 

OUTPUT 

The 
that  is 
entered 
defined 

remaining  input  data  to  the  diesel  program  are  contained 
read,  format  free,  from  an  input  file  titled  DSLINP.  The 
as  shown  in  table  I,  where  the  variables  in  the  input  data 
as  follows: 

in  a 
data 
set 

data 

must 

are 

set 

be 

RPM 

diesel  engine  crankshaft  speed,  rev/min 

WFCY 

weight  (mass)  of  fuel  added  per  cycle,  lbm 

BORE 

cylinder  bore,  in. 

STROKE 

stroke,  in. 

CON ROD 

connecting  rod  length,  in. 

CR 

geometric  compression  ratio 

MODE 

piston  cylinder  mode  (MODE  -  1  for  conventional  piston-cylinder 
configuration,  MODE  =  2  for  opposed  piston  configuration) 

Cl 

shape  parameter  for  diffusion  burning 

C2 

crank  angle  for  start  of  combustion 

C3 

combustion  duration  parameter 

C4 

premixed  burning  fraction 

C5 

shape  parameter  for  premixed  burning 

ANNND 

multiplier  for  Annand  heat  transfer  correlation 

THEAD 

cylinder  head  temperature,  °R 

TPISTN 

piston  temperature,  °R 

TSLEEV 

cylinder  wall  temperature,  °R 

MEXH 

exhaust  valve  or  port  indicator  (MEXH  a  0  for  exhaust 
MEXH  =  1  for  exhaust  valves) 

ports , 

EVO 

exhaust  valve  or  port  opening  crank  angle,  °ATDC 

EVC 

exhaust  valve  or  port  closing  crank  angle,  °ATDC 
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MINT 

AVO 

AVC 

NTEXH 

CDEXH 

ALPHEX(I) 

FEXH(I) 

WIDTHE 

NTINT 

CDINT 

ALPHIN(I) 

FINT(I) 

WIDTHI 


intake  valve  or  port  indicator  (MINT  =  0  for  intake  ports,  MINT  =  1 
for  Intake  valves) 

intake  valve  or  port  opening  crank  angle,  °ATDC 

Intake  valve  or  port  closing  crank  angle,  °ATDC 

number  of  exhaust  valve  area  versus  crank  angle  data  points  to  be 
read  in 

exhaust  valve  or  port  flow  coefficient 

array  of  crank  angles  for  exhaust  valve  flow  areas,  °ATDC 

array  of  exhaust  valve  flow  areas,  in.2 

fraction  of  the  cylinder  circumference  devoted  to  exhaust  port 

number  of  intake  valve  area  versus  crank  angle  data  points  to  be 
read  in 

intake  valve  or  port  flow  coefficient 

array  of  crank  angles  for  Intake  valve  flow  areas,  °ATDC 

array  of  intake  valve  flow  areas,  in.2 

fraction  of  the  cylinder  circumference  devoted  to  intake  port 


When  MEXH  =  1,  exhaust  valves  are  being  used  and  the  user  must  specify 
NTEXH,  CDEXH,  and  the  arrays  of  crank  angle  versus  flow  area.  When  MEXH  =  0, 
exhaust  ports  are  used  and  the  user  need  only  specify  WIDTHE  and  CDEXH.  The 
case  Illustrated  previously  corresponds  to  an  engine  that  has  exhaust  valves 
and  intake  ports. 


OUTPUT 

The  output  from  the  program  simulating  a  two-stroke  diesel  engine  consists 
of  two  parts.  The  first  part,  a  summary  of  calculated  quantities  that  charac¬ 
terize  the  engine's  performance,  is  discussed  in  detail  hereinafter.  The  sec¬ 
ond  part  is  a  listing  of  the  cylinder  pressure  and  temperature,  and  other 
quantities  at  5°  Intervals  during  the  engine  cycle.  The  abbreviations  used  to 
identify  the  columns  are  defined  as  follows: 

CA  crank  angle  at  which  data  were  calculated;  top  dead  center  corre¬ 

sponds  to  0° 

PCYL  cylinder  pressure 
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TCYL  cylinder  temperature.  Since  the  program  is  based  on  a  single-zone 
combustion  model,  a  single  temperature  is  defined  for  the  entire 
cyl inder  contents. 

PHICYL  equivalence  ratio  of  the  cylinder  gases.  This  is  the  ratio  of  the 
fuel-to-air  mass  ratio  in  the  cylinder  to  the  stoichiometric  fuel- 
to-air  mass  ratio. 

MCYL  mass  of  cylinder  contents 

MDOTIN  mass  flow  rate  into  cylinder 

MDOTOUT  mass  flow  rate  leaving  the  cylinder 

BALANCE  the  amount  of  imbalance  in  the  energy  balance.  It  is  calculated  and 
listed  as  a  check  on  the  solution  of  the  differential  equations. 

If  changes  are  made  to  the  program,  such  as  the  addition  of  new 
subroutine  modules,  the  energy  balance  should  be  monitored  to 
ensure  that  it  stays  close  to  zero. 

The  summary  part  Includes  a  listing  of  the  input  conditions  for  the  run 
and  gives  the  engine  configuration  and  the  engine  operating  condition.  Then  a 
listing  of  calculated  quantities  is  provided  based  on  the  program  results.  The 
first  quantities  calculated  are  the  air  and  fuel  flow  rates  per  cylinder.  The 
flow  rates  are  calculated  by  dividing  the  total  amount  of  air  and  fuel  added  to 
the  engine  over  one  cycle  by  the  time  required  for  one  cycle. 

The  IMEP  is  calculated  by  summing  the  work  done  over  the  cycle  and  then 
dividing  by  the  displacement  volume.  The  BMEP  is  then  calculated  by  subtract¬ 
ing  the  FMEP  that  was  calculated  by  using  an  empirical  correlation  from  the 
IMEP.  The  same  process  is  used  to  obtain  the  indicated  and  brake  power,  the 
indicated  and  brake  specific  fuel  consumption,  and  the  Indicated  and  brake 
thermal  efficiency. 

During  the  integration  process  of  solving  the  differential  equations  of 
the  mathematical  model,  the  maximum  values  of  cylinder  temperature  and  pressure 
are  recorded  along  with  the  time  at  which  they  occur.  These  results  are  also 
provided  In  the  summary  section. 

In  the  literature  a  variety  of  parameters  have  been  defined  to  character¬ 
ize  the  scavenging  process.  These  quantities  can  be  calculated  from  a  know¬ 
ledge  of  how  much  residual  gas  Is  present  In  the  cylinder  and  how  much  fresh 
charge  was  added.  Since  the  Instantaneous  equivalence  ratio  In  the  cylinder  Is 
known  from  the  solution  of  equation  <6),  we  can  determine  the  values  of  $  at 
the  point  where  the  exhaust  valve  opens,  $evo»  anc*  at  the  P<^nt  where  compres¬ 
sion  starts,  ♦cmp.  The  total  mass  In  the  cylinder  Is  also  known  from  solution 
of  the  cylinder  mass  balance.  To  characterize  the  cylinder  contents,  It  Is 
convenient  to  divide  the  total  cylinder  mass  Into  the  following  five 
categories: 
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mFC 

mass 

of 

mfuel ,new 

mass 

of 

mR,fuel 

mass 

of 

mR,BA 

mass 

of 

mR,UA 

mass 

of 

fresh  charge  (pure  air) 
new  fuel  injected 
burned  fuel  In  residual 
burned  air  In  residual 
unburned  air  in  residual 


The  definitions  of  these  quantities  are  based  on  the  assumption  that  the 
equivalence  ratio  in  the  cylinder  is  less  than  1.0.  This  is  a  good  assumption 
for  diesel  engines  since  they  normally  operate  with  a  maximum  equivalence  ratio 
of  0.6  to  0.8.  The  definitions  also  assume  that  the  combustion  reactions  cause 
the  fuel  and  air  mixture  to  go  to  equilibrium  products.  These  quantities  can 
be  determined  by  solving  the  following  five  equations: 


*  f  = - - - 

cmp  s  mFC  +  mR  BA  +  mR  yA 

.  p  mR,fuel  *  mfuel .new 
evo  s  mFC  +  mR>BA  +  mR  m 

Vfuel  =  F 
mR,BA  S 

A  F  „  mR.fue1 
evo  s  mR,BA  +  mR,UA 

mR,fuel  +  mR,BA  +  mR,UA  +  mFC  +  mfuel,new  =  mcyl 
These  equations  can  be  reduced  to 


(13a) 

(13b) 

(13c) 

(13d) 

(13e) 


cmp 

mR.BA  *  l  +  <t»evoFs  mcyl 

(14a) 

mR,fuel  *  FsmR,BA 

(14b) 

mR,UA  "  (1/*evo  "  1)mR,BA 

(14c) 

mFC  ’  <1/4>cmp  "  ^  /^evo)mR,BA 

( 14d) 

'fuel, new  ”  ^evo/(^cmp  “  1)FsmR,BA 

(14e) 
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Once  the  mass  In  the  cylinder  has  been  charactered,  the  various  scavenging 
Indicators  can  be  calculated. 


purity  =  P  = 


mFC  *  VUA 


mFC  +  mR,BA  +  mR,UA  +  mR,fuel 


mR  BA  +  mR  tlA  +  mR  fuel 

residual  fraction  =  f  =  - - ^ ^ - *4^ - 

mFC  +  mR,BA  +  mR,UA  +  mR,fuel 

mFC 

trapping  efficiency  =  n^r  =  total  mass  of  air  entering  engine 

mr 


charging  efficiency  =  nch 


TC 


p1ntakevd1 sp 


(15) 

(16) 

(17) 

(18) 


del  1  very  ratio  = 


(19) 


The  scavenging  efficiency,  which  Is  defined  as  the  mass  of  delivered  air 
that  Is  retained,  divided  by  the  mass  of  trapped  cylinder  charge.  Is  calculated 
using  the  following  equation  that  Is  developed  In  Schweitzer  (ref.  4): 


scavenging  efficiency  =  n 

j 


1 


1  + 


(1/P  -  1) 


4>, 


evo 


(20) 


The  summary  also  determines  what  happened  to  the  fuel  energy.  The  frac¬ 
tion  of  the  fuel  energy,  as  reflected  by  the  lower  heating  value,  which  goes 
to  work,  to  heat  transfer,  and  to  the  exhaust.  Is  calculated 


1  work  = 


’tot 


mfuel ,newHVL 


%  heat  loss 


'tot 


mfuel ,newHVL 


1  exhaust  =  1  -  X  work  -  1  heat  loss 


(21) 

(22) 

(23) 


where  W+0t  Is  the  total  amount  of  work  done  during  the  cycle,  and  Qtot 
the  total  heat  transfer  from  the  cylinder  gases  to  the  wall,  and  HVL  Is  the 
lower  heating  value  of  the  fuel. 
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Finally,  the  summary  provides  an  estimate  of  the  engine's  exhaust  tempera¬ 
ture.  Although  the  actual  properties  of  the  mass  leaving  the  exhaust  port  are 
time  varying,  the  value  of  temperature  that  provides  the  same  average  enthalpy 
transfer  as  the  time-varying  flow  Is  determined  by  solving  the  following  equa¬ 
tion  Iteratively: 


^exh’^exh’^exh^ 


J  houtdmout 

J  d"W 


where 

Texh  exhaust  temperature  (unknown) 

Pexh  exhaust  manifold  pressure 

<j>exh  equivalence  ratio  In  exhaust 

This  provides  a  mass-averaged  exhaust  temperature. 


(24) 


SAMPLF  CASE 

As  a  sample  case,  consider  a  two-stroke  diesel  engine  with  the  following 
characteristics: 

opposed  piston  design 

0.00022  lbm  fuel  Injected  per  stroke 

150  psla  Intake  pressure 

885  °R  Intake  air  temperature 

137.5  psla  exhaust  pressure 

3.101  In.  bore  by  2.940  In.  stroke 

7.168  In.  connecting  rod  length 

9.1712  compression  ratio 

Intake  port  opens  at  125  °ATDC,  occupies  55%  of  cylinder  bore,  flow 
coefficient  of  0.8 

exhaust  port  opens  at  100  °ATDC,  occupies  55%  of  cylinder  bore,  flow 
coefficient  of  0.8 

The  program  below  calls  the  subroutine  DSL2  to  calculate  the  performance 
of  the  two-stroke  diesel  engine.  The  Intake  pressure  and  temperature  and  the 
exhaust  pressure  are  Input  through  the  subroutine  argument  list.  The  remaining 
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data  are  provided  In  the  Input  data  set  named  DSL.INP.  This  data  set  Is  provi¬ 
ded  in  table  II.  The  complete  output  file  Is  also  printed  below. 

C  SAMPLE  PROGRAM 

C 

PI-150. 

T 1=885 . 

PEXH=1 37 . 5 

C 

CALL  DSL2(TI ,PI ,PEXH,TEXH,XM) 

C 

STOP 

END 


Complete  Output  File 


BORE 

STROKE 

CONNECTING  ROD 

GEOMETRIC  COMPRESSION  RATIO 

EFFECTIVE  COMPRESSION  RATIO 

SWEPT  VOLUME 

ENGINE  SPEED 

MEAN  PISTON  SPEED 

SUPPLY  AIR  PRESSURE 

SUPPLY  AIR  TEMPERATURE 

EXHAUST  PRESSURE 


3.101  INCHES 
2.940  INCHES 
7.168  INCHES 
9.17 
6.21 

44.409  CUBIC  INCHES 
6122.  RPM 
3000.  FT/M IN 
150.00  PSIA 
885.00  DEG  R 
137.50  PSIA 


PORT  TIMING:  INTAKE  OPEN  125. 

INTAKE  CLOSE  235. 
EXHAUST  OPEN  100. 
EXHAUST  CLOSE  260. 


FUEL  FLOW  RATE 

0.0224  LBM/SEC 

FUEL  INJECTED/CYCLE 

0.000220 

LBM/CYCLE 

AIR  FLOW  RATE 

0.6544  LBM/SEC 

AIR  INDUCTED/CYCLE 

0.006413 

LBM/CYCLE 

IMEP 

314.7  PSI 

IHP 

216.1  HP 

BMEP 

286.7  PSI 

BHP 

196.9  HP 

FMEP 

28.0  PSI 

FHP 

19.2  HP 

ISFC  0.3739  LBM/HP-HR 

BSFC  0.4104  LBM/HP-HR 

INDICATED  THERMAL  EFFICIENCY  0.3738 

BRAKE  THERMAL  EFFICIENCY  0.3405 
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MAXIMUM  PRESSURE 
MAXIMUM  TEMPERATURE 
MAXIMUM  RATE  OF  PRESS  RISE 


PURITY  0.7343 
SCAVENGING  EFFICIENCY  0.6965 
TRAPPING  EFFICIENCY  0.5941 
CHARGING  EFFICIENCY  0.3241 
DELIVERY  RATIO  0.5455 
RESIDUAL  FRACTION  0.3164 


2784.4  PSIA  AT  368.8  DEG 

4571.5  DEG  R  AT  378.0  DEG 

82.81  PSI/DEG  AT  355.0  DEG 


EQUIVALENCE  RATIO  DURING  COMPRESSION  0.2528 
EQUIVALENCE  RATIO  AT  EVO  0.8304 
EQUIVALENCE  RATIO  BASED  ON  FUEL  AND 

AIR  FLOW  RATES  (EXHAUST)  0.4932 


% 

FUEL 

ENERGY 

TO 

HEAT 

LOSS 

8.11 

« 

% 

FUEL 

ENERGY 

TO 

WORK 

37.38 

% 

% 

FUEL 

ENERGY 

TO 

EXHAUST 

54.52 

% 

EXHAUST 

1  TEMPERATURE 

2057 

.5  DEG 

R 

CA 

PCYL 

TCYL 

PHICYL 

MCYL 

MDOTIN 

MDOTOUT 

105. 

313.0 

3068.1 

0.8304 

0.005692 

0.000000 

0.000040 

110. 

276.0 

2985.0 

0.8304 

0.005409 

0.000000 

0.000071 

115. 

236.2 

2885.8 

0.8304 

0.005000 

0.000000 

0.000090 

120. 

197.9 

2776.6 

0.8304 

0.004528 

0.000000 

0.000096 

125. 

165.3 

2669.0 

0.8304 

0.004075 

0.000000 

0.000082 

130. 

142.8 

2578.1 

0.8271 

0.003762 

0.000012 

0.000043 

135. 

137.4 

2512.4 

0.8014 

0.003818 

0.000030 

-0.000008 

140. 

137.5 

2446.6 

0.7624 

0.004022 

0.000043 

0.000000 

145. 

137.6 

2366.5 

0.7162 

0.004249 

0.000055 

0.000007 

150. 

137.7 

2279.5 

0.6666 

0.004496 

0.000065 

0.000014 

155. 

137.9 

2190.6 

0.6168 

0.004756 

0.000073 

0.000020 

160. 

138.1 

2103.2 

0.5689 

0.005021 

0.000079 

0.000026 

165. 

138.3 

2019.8 

0.5241 

0.005286 

0.000084 

0.000032 

170. 

138.6 

1941.8 

0.4830 

0.005542 

0.000087 

0.000037 

175. 

138.8 

1870.2 

0.4459 

0.005785 

0.000089 

0.000042 

180. 

139.0 

1805.2 

0.4127 

0.006010 

0.000089 

0.000046 

185. 

139.2 

1746.8 

0.3833 

0.006210 

0.000087 

0.000050 

190. 

139.4 

1694.7 

0.3575 

0.006383 

0.000084 

0.000053 

195. 

139.6 

1648.8 

0.3350 

0.006525 

0.000080 

0.000055 

200. 

139.8 

1608.7 

0.3156 

0.006632 

0.000074 

0.000056 

205. 

140.0 

1574.3 

0.2990 

0.006702 

0.000067 

0.000057 

210. 

140.1 

1545.1 

0.2851 

0.006733 

0.000059 

0.000056 

215. 

140.3 

1521.2 

0.2737 

0.006723 

0.000049 

0.000055 

220. 

140.5 

1502.4 

0.2647 

0.006672 

0.000038 

0.000053 

225. 

140.7 

1488.8 

0.2582 

0.006578 

0.000027 

0.000050 

230. 

140.9 

1480.5 

0.2542 

0.006441 

0.000014 

0.000045 

235. 

141.1 

1478.0 

0.2528 

0.006261 

0.000000 

0.000040 

240. 

141.9 

1480.1 

0.2528 

0.006070 

0.000000 

0.000037 

245. 

143.8 

1485.0 

0.2528 

0.005894 

0.000000 

0.000033 

250. 

147.0 

1493.3 

0.2528 

0.005739 

0.000000 

0.000028 

BALANCE 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 
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255.  152.4  1506.9  0.2528  0.005622  0.000000  0.000018  0.00000 
260.  161.5  1529.1  0.2528  0.005574  0.000000  0.000000  0.00000 
265.  174.3  1558.4  0.2528  0.005574  0.000000  0.000000  0.00000 
270.  189.3  1590.8  0.2528  0.005574  0.000000  0.000000  0.00000 
275.  207.1  1626.5  0.2528  0.005574  0.000000  0.000000  0.00000 
280.  228.1  1665.9  0.2528  0.005574  0.000000  0.000000  0  00000 
285.  253.3  1709.2  0.2528  0.005574  0.000000  0.000000  0 1 00000 
290.  283.5  1756.8  0.2528  0.005574  0.000000  0.000000  0.00000 
295.  320.0  1809.0  0.2528  0.005574  0.000000  0.000000  0^00000 
300.  364.2  1866.2  0.2528  0.005574  0.000000  0.000000  0.00000 
305.  417.9  1928.8  0.2528  0.005574  0.000000  0.000000  0.00000 
310.  483.6  1996.9  0.2528  0.005574  0.000000  0.000000  0.00000 
315.  563.9  2070.7  0.2528  0.005574  0.000000  0.000000  0.00000 
320.  661.7  2149.8  0.2528  0.005574  0.000000  0.000000  0.00000 
325.  780.0  2233.5  0.2528  0.005574  0.000000  0.000000  0.00000 
330.  920.6  2320.5  0.2528  0.005574  0.000000  0.000000  0.00000 
335.  1082.9  2408.2  0.2528  0.005574  0.000000  0.000000  0.00000 
340.  1261.4  2492.7  0.2528  0.005574  0.000000  0.000000  0.00000 
345.  1442.5  2568.6  0.2528  0.005574  0.000000  0.000000  0.00000 
350.  1704.3  2788.6  0.2848  0.005586  0.000000  0.000000  0.00000 
355.  2094.0  3230.4  0.3707  0.005619  0.000000  0.000000  0.00000 
360.  2481.3  3730.8  0.4848  0.005662  0.000000  0.000000  0.00000 
365.  2730.4  4151.8  0.5986  0.005706  0.000000  0.000000  0.00000 
370.  2776.7  4426.4  0.6916  0.005741  0.000000  0.000000  0.00000 
375.  2641.3  4553.0  0.7563  0.005766  0.000000  0.000000  0.00000 
380.  2392.7  4565.5  0.7950  0.005781  0.000000  0.000000  0.00000 
385.  2100.3  4504.7  0.8152  0.005788  0.000000  0.000000  0.00000 
390.  1811.6  4403.1  0.8245  0.005792  0.000000  0.000000  0.00000 
395.  1550.9  4282.9  0.8283  0.005793  0.000000  0.000000  0.00000 
400.  1326.8  4157.7  0.8296  0.005794  0.000000  0.000000  0.00000 
405.  1138.8  4034.6  0.8300  0.005794  0.000000  0.000000  0.00000 
410.  983.1  3917.2  0.8302  0.005794  0.000000  0.000000  0.00000 
415.  854.8  3807.0  0.8302  0.005794  0.000000  0.000000  0.00000 
420.  748.9  3704.4  0.8302  0.005794  0.000000  0.000000  0.00000 
425.  661.5  3609.6  0.8302  0.005794  0.000000  0.000000  0.00000 
430.  588.9  3522.2  0.8302  0.005794  0.000000  0.000000  0.00000 
435.  528.5  3441.8  0.8302  0.005794  0.000000  0.000000  0.00000 
440.  477.8  3368.1  0.8302  0.005794  0.000000  0.000000  0.00000 
445.  435.2  3300.5  0.8302  0.005794  0.000000  0.000000  0.00000 
450.  399.2  3238.7  0.8302  0.005794  0.000000  0.000000  0.00000 
455.  368.6  3182.3  0.8302  0.005794  0.000000  0.000000  0.00000 
460.  342.6  3130.7  0.8302  0.005794  0.000000  0.000000  0.00000 
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CONCLUSION 


The  program  described  in  this  report  can  be  used  to  simulate  any  two- 
stroke  engine  with  either  a  standard  piston-cylinder  configuration  or  an 
opposed  piston  configuration.  The  program  was  developed  to  provide  a  way  to 
predict  the  performance  of  the  diesel  core  of  a  high-output  combined-cycle 
diesel  engine.  The  sample  case  presented  is  one  example  of  a  possible  engine 
suitable  for  this  use. 

This  program  can  be  used  to  investigate  the  effects  on  engine  performance 
of  combustion  timing,  valve  or  port  timing,  heat  transfer,  and  engine  geometry 
Gross  engine  parameters  such  as  bore,  stroke,  connecting  rod  length,  and  com¬ 
pression  ratio  as  well  as  engine  operating  parameters  such  as  speed,  load,  and 
altitude  can  be  easily  varied.  While  exact  prediction  of  engine  characteris¬ 
tics  is  not  possible,  the  program  should  accurately  predict  trends. 
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TABLE  I.  -  INPUT  VARIABLES 


TABLE  II.  -  DSL. INP 


6122. 

0.00022 

3.101 

2.940 

7.168 

9.1712 

2 

1 .0 

345. 

55. 

0.0 

3.5 

0.7 

1460. 

1460. 

1460. 

0 

100. 

260. 

0 

125. 

235. 

0.55 

0.8 

0.55 

0.8 
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